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Abstract 

We compute the real part of the 2-loop vertex corrections for charmless hadronic 
B decays, completing the NNLO calculation of the topological tree amplitudes in 
QCD factorization. Among the technical aspects we show that the hard-scattering 
kernels are free of soft and collinear infrared divergences at the 2-loop level, which 
follows after an intricate subtraction procedure involving evanescent four quark 
operators. The numerical impact of the considered corrections is found to be mod- 
erate, whereas the factorization scale dependence of the topological tree amplitudes 
is significantly reduced at NNLO. We in particular do not find an enhancement of 
the phenomenologically important ratio \C/T\ from the perturbative calculation. 
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1 Introduction 



The study of hadronic B meson decays into a pair of light (charmless) mesons reveals 
interesting information about the underlying four quark interactions and the related phe- 
nomenon of CP violation. While these decay modes are intensively investigated at current 
and future B physics experiments, the main challenge for precise theoretical predictions 
consists in the computation of the hadronic matrix elements. QCD factorization [I], or 
its field theoretical formulation in the language of Soft-Collinear Effective Theory [2], 
is a systematic framework to compute these matrix elements from first principles. The 
starting point is a factorization formula, which holds in the heavy quark limit m& — > oo, 

(M 1 M 2 \Q i \B) ~ F* Ml (0)f Ma J du T^u)<p M2 {u) 

+ Ib fhh fhi 2 J dudvdu T t n (u, v, u) <f> B (w) 4>Ah 0) 4>m 2 (u), (1) 

where the perturbatively calculable hard-scattering kernels T?' 11 encode the short- distance 
strong-interaction effects and the non-perturbative physics is confined to some process- 
independent hadronic parameters such as decay constants /m, light-cone distribution 
amplitudes <$m and a transition form factor F+ Ml at maximum recoil q 2 = 0. 

In this work we address perturbative corrections to the factorization formula ([1]). 
Whereas next-to-leading order (NLO) corrections to the hard-scattering kernels T*' 11 are 
known from the pioneering work in pQ, partial next-to- next-to-leading order (NNLO) 
corrections have recently been worked out [31 HI El E]- The a 2 s corrections to the kernels 
T/ {spectator scattering) are by now completely determined to NNLO: the corrections for 
the topological tree amplitudes have been computed in [3] and the ones for the so-called 
penguin amplitudes in [3]. 

In contrast to this the computation of o? s corrections to the kernels T/ (vertex cor- 
rections) is to date incomplete. Whereas we computed the imaginary parts of the hard- 
scattering kernels for the topological tree amplitudes in [5J E], we complete the NNLO 
calculation of the tree amplitudes in this work by computing the respective real parts. 
Partial results of this calculation, in particular the analytical expressions of the required 
2-loop Master Integrals, have already been given in [6]. 

The organization of this paper is as follows: The technical aspects of the NNLO 
calculation are presented in Section [2J We start by briefly recalling our definitions and 
conventions and make some remarks concerning the computation of the 2-loop diagrams. 
We then show in some detail how to extract the hard-scattering kernels from the matrix 
elements which are formally infrared divergent. This subtraction procedure, which be- 
comes particularly involved for the colour-suppressed tree amplitude, is complicated due 
to the presence of evanescent four quark operators which arise in intermediate steps of 
the calculation. Our analytical results for the hard-scattering kernels are summarized in 
Section El We briefly discuss the numerical impact of the considered NNLO corrections 
in Section HI before we conclude in Section El Several technical issues of the calculation 
and the explicit expressions of the hard-scattering kernels are relegated to the Appendix. 
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2 NNLO calculation 



The calculation of the real parts of the topological tree amplitudes proceeds along the 
same lines as the one of the imaginary parts that we presented in [5]. Still, the current 
calculation turns out to be considerably more complex in several respects. First, it 
requires the calculation of a larger amount of 2-loop integrals, which are in addition 
more complicated since they involve up to three (instead of one) massive propagators. 
Second, the renormalization procedure and the infrared subtractions reveal their full 2- 
loop complexity only in the current calculation as a consequence of the fact that the 
tree level contribution is real. In the following we summarize the technical aspects of 
the calculation and refer for a more detailed description of the general strategy to [5] 
(cf. also [6]). 



2.1 Operator basis 

The topological tree amplitudes can be derived from the hadronic matrix elements of the 
current-current operators in the effective weak Hamiltonian 

H cS = ^ V* d V ub (CxQi + C 2 Q 2 ) + h.c. (2) 

As we apply Dimensional Regularization3 (DR) to regularize ultraviolet (UV) and infrared 
(IR) singularities, evanescent four-quark operators appear in intermediate steps of the 
calculation. The full operator basis required for the present calculation becomes^ 

Q x = [u^LT A b] [d lp LT A u], 
Q 2 = [uYLb] [d^Lu], 

E x = [uYl v l p LT A b] [d llllulp LT A u] - 16 Q 1; 
E 2 = [u^YYLb] [dwypLu] - 16 Q 2 , 

E[ = [uYYYYYL T A b] [d lp ri ulplalr L T A u] -20 Ex- 256 Q h 

E' 2 = [fi/yy/fi b] [d llllvlplalr Lu] —20E 2 — 256 Q 2 , (3) 

with colour matrices T and L = 1 — 75. We stress that previous studies within QCD 
factorization, as e.g. [U El Hj, have often been formulated in a different operator basis 
with a Fierz-symmetric definition of the physical operators. As has been argued in [5], it 
is more convenient for the current calculation to use the operator basis Q since it allows 
to work with a naive anticommuting 75 beyond NLO [7]. 

There are two different insertions of a four-quark operator which are illustrated in 
Figure 1 of [5j. The first one gives rise to the colour- allowed tree amplitude ai(MiM 2 ), 
which corresponds to the flavour content [q s b] of the decaying B meson, [q s u] of the recoil 



1 We write d — 4 — 2e and use an anticommuting 75 according to the NDR scheme. 
2 This operator basis has been named CMM basis in 5 (denoted by a hat). 
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meson Mi and [ud] of the emitted meson M 2 . The colour-suppressed tree amplitude 
a2{MiM 2 ) follows from the second insertion and belongs to the flavour contents [q s b], 
[q s d] and [««], respectively. In [5] we did not consider the second type of insertions since 
we could derive the imaginary part of the colour-suppressed amplitude from the one of 
the colour-allowed amplitude using Fierz-symmetry arguments^. 

In the current calculation we cannot proceed along the same lines, since a Fierz- 
symmetric operator basis has not yet been worked out to NNLO0. We therefore consider 
both types of insertions in this work, which also provides an independent cross-check of 
our previous result for the imaginary part of the colour-suppressed tree amplitude. 



2.2 2-loop calculation 

The main task of the calculation consists in the computation of a large number of 2-loop 
diagrams (shown in Figure 2 of [5]). We use an automatized reduction algorithm, which 
is based on integration- by-parts techniques [11] , to express these diagrams in terms of an 
irreducible set of Master Integrals (Mis). In addition to the Mis that appeared in the 
calculation of the imaginary part of the NNLO vertex corrections (cf. Figure 3 of [5]), we 
find 22 Mis which are shown in Figure [H In total the current calculation requires the 
computation of 36 Mis to up to five orders in the e-expansion. 

Apart from the Mis that involve the charm quark mass, the analytical results for the 
Mis from Figure [T] can be found in [BJT . The Mis can be expressed in terms of Harmonic 
Polylogarithms (HPLs) [14] of weight w < 4, 



H(0;x) 


= ln(*), 


#(0,0,1;*) 


= Li 3 (*), 


H(l;x) 


= — ln(l — x), 


#(0,1,1;*) 


= Si, 2 (*) 


H(-l;x) 


= ln(l + x), 


#(0,0, 0,1; a;) 


= Li 4 (x), 


#(0,1; x) 


= Li 2 (x), 


#(0,0,1,1;*) 


= S 2 , 2 (*) 


#(0,-1;*) 


= -Li 2 (-*), 


#(0,1,1,1;*) 


= Si,s(aO 


#(-1,0,1;*) 


= Wi(z), 


#(0,-1,0,1;*) 


= n 2 (x), 



(4) 



where we introduced a shorthand notation for the last two HPLgj. Moreover, the massive 
non-planar 6-topology MI (last diagram from Figure [1]) involves a constant in the finite 
term which, until recently, was only known numerically, Cq = —60.2493267(10) [13] . In a 
recent work it was shown that its analytical value is Co = — 1677r 4 /270 [15] . 

3 To do so we introduced a second operator basis named traditional basis in |5J (denoted by a tilde). 
4 We emphasize that the operator basis from Section 8 in [5] is not Fierz-symmetric and the one from 
Appendix A in [9] is presumably not either [TO] , 

5 Part of these results have recently been confirmed by various groups [T2l [13] . 

6 The explicit expression of 7ii(x) in terms of Nielsen Polylogarithms can be found e.g. in equation 
(10) of [15] . On the other hand Tiiix) has to be evaluated numerically (in Section [3?2l we find, however, 
analytical expressions in the convolutions with the light-cone distribution amplitude of the meson M^). 
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Figure 1: Additional Master Integrals that appear in the calculation of the real 
parts of the NNLO vertex corrections. Dashed/double/wavy internal lines denote 
propagators with mass / / m c . Dashed/solid/double external lines correspond to 
virtualities / um\ / m\ . Dotted propagators are taken to be squared. 

The charm mass dependent Mis can be found in [6l [15]. In this case there exist 
analytical results apart from the finite terms of two 4-topology Mis. We may, however, 
evaluate these contributions numerically to implement charm mass effects in the current 
analysis. 

2.3 Renormalization 

The calculation of the renormalized matrix elements requires standard counterterms from 
QCD and the effective Hamiltonian. We write the renormalized matrix elements as 

(Qi) = %ij (Qj)bare; (5) 

where contains the wave-function renormalization factors of the quark fields and Z is 
the operator renormalization matrix in the effective theory. Here and below we introduce 
a shorthand notation for the perturbative expansions, 

00 k °° k 

(Qi)(bare) = £ (£) <Q«>gL>> % = ** + £ (S) ^ " ^ 

k=0 k=l 

It turns out that the wave-function renormalization factors in can be neglected in the 
calculation of the hard-scattering kernels since they are absorbed by the form factor and 
the light-cone distribution amplitude in the factorization formula, which are defined in 
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terms of full QCD fields (rather than HQET or SCET fields), for details cf. Section 4.2 
of [5]. We renormalize the coupling constant in the MS-scheme, 



e 7B /i 2 \ £ „, , 3 -2e 



and the 6-quark mass in the on-shell scheme, 
The 1-loop and 2-loop MS operator renormalization matrices can be inferred from [8j [17] 



2 i A 2 n \ 1 

Z 3 12 g U U \ 1 



Z (2) 



6 1 J e 

17_2„ _26 , 4 _25 , _5_„ _31 i J_„ 19 5 \ 

1( 3' 4 / 3 " r 9 "7 6 " r 36'V 18" r 27' t / 96 108 I ±_ 

-39 + 271/ 4 -f + £», A J J 

79 , 4 _205,10 1531 _J^ n „ _±_± Uf JL _35_ > , 

12 ' 9 J 18 ' 27 J 288 216 J 72 81 / 384 864 1 1 



_|_ I 12 1 9 J 18 1 27 J 288 216 J 72 81 J 384 864 _ /g\ 

1 83 + 5 o 119_J_„, 8 _^5_ _J_ J e ' 

4^3"/ d 16 18 "7 9 192 72 



where the lines refer to the physical operators and the columns to the full operator basis 
including the evanescent operators from ([3]). 

2.4 IR subtractions 

In order to extract the hard-scattering kernels Tj we rewrite the renormalized matrix 
elements in the factorized form 

(Qi) = F ■ Ti ® $ + . . . (10) 

where F denotes the form factor, $ the product of decay constant and distribution 
amplitude, (g> the convolution integral and the ellipsis the spectator scattering term which 
we disregard in the following. As has been discussed in detail in Section 4.2 of [5], only 
naively non-factorizable (nf) 1-loop diagrams contribute to the NLO kernels, 

(Q^ + Z^( Qj }^ = F^.T^®^. (11) 

Similarly, the calculation of the NNLO kernels involves only non-factorizable 2-loop dia- 
grams (but factorizable (f) 1-loop diagrams), 

wS , +4 1) [(^■)S ) + (^ 1) ]+^ ) (Qi) (0) 

= F<°> • 7f ® $(°> + F« p • ^ ® + F<°) • lf> ® $« p , (12) 
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where the subscript " amp" (amputated) has been introduced to denote corrections with- 
out wave-function renormalization. We see that the calculation of the NNLO kernels re- 
quires the NLO kernels to 0(e 2 ) as they enter (fT2]) in combination with the IR-divergent 
form factor correction F&mp ~ 1/ £ ir- As a consequence the factorization formula has 
to be extended in intermediate steps of the calculation to include evanescent operators, 
which have to be renormalized such that their (IR-fmite) matrix elements vanish (for 
details cf. Section 4.3 of [5]). 

At NNLO the subtraction procedure becomes somewhat involved. It is particularly 
complicated in the calculation of the colour-suppressed tree amplitude, where a Fierz- 
evanescent operator appears at tree level. In the following we discuss the subtraction 
procedure in some detail. Throughout this section we concentrate on the real parts of 
the hard-scattering kernels, since the respective imaginary parts have already been given 
in [5]. We refer to Appendix |X] for the explicit expressions of the auxiliary coefficient 
functions ti(u) that we introduce below. 

Colour-allowed tree amplitude 

To NNLO we find three operators that contribute to the right hand side of fllOp . In 
the position space representation they correspond to products of a local heavy-to-light 
current u(x)Tib(x) and a non-local light-quark current d(y)[y,x]T 2 u(x), where the usual 
gauge link factor [y, x] is understood. We choose the basis of Dirac structures Ti ®T 2 as3 

O = [fL] ® [ 7/l L] , 

o E = [YYYL] ® [i^iuiA - 16 o, 

O e > = [YYYYYL] ® [ 7At 7.7p7<x7r L] -20O E - 256 O, (13) 

such that the factorized hadronic matrix element of O gives the standard QCD form 
factor and the light-cone distribution amplitude of the emitted meson M 2 . The operators 
Oe and Oe> are evanescent. 

We first compute (ITT]) to 0(e 2 ) to determine the NLO kernels. We find that the 
colour-singlet kernels vanish, T 2 = T^e = T 2 e> = 0> while the colour-octet kernels 
become 

Re Ti x \u) = ^ (^j |to(«) - + (t^u) + 3L 2 )e + (t 2 (u) - L')e 2 + 0(e 3 )}, 
Re T« (u) = (j^y {^,oH + 2L+ (t E ^u) - L 2 )e + 0(e 2 )}, (14) 

and T^ E , = with L = ln/i 2 /m 2 . The IR subtractions on the right hand side of (fT2l 
require in addition form factor and wave function corrections to the operators O and O e 
(they can be found in Section 4.3 of [5]). We finally perform the convolutions of the NLO 

7 We do not consider colour-octet operators since their hadronic matrix elements vanish. 
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kernels with the wave function corrections, which yields 

^ (0) Re If $« p = f c (Q £ {^ + UM + 0(e)} F® $C°) (15) 
and an additional //-dependent contribution to the physical kernel from 

F<°> Re T« $« P> * - 3 |i2L + t E » + 0(e)} F® $(°). (16) 

Colour-suppressed tree amplitude 

In this case we find an analogous set of operators, 

6 = [<fL] ® [i,L] , 
O e = [YYYL] ® [ 7 ^ 7l , 7p L] - 16 0, 

B , = [ 7 ^ 7 - 7 ^ 7 - 7 -L] (§) [7 M 7,7p7<x7. L] -20O E - 256 0, (17) 

but the fields are now given in the wrong ordering u(y)[y,x]Tib(x) and d(x)T2u(x) (indi- 
cated by (§)), which does not yield a form factor and a light-cone distribution amplitude. 
The latter follow from the factorized hadronic matrix element of the operator 

O = d(x)YLb(x) <g> u(y)[y,x]-ff,Lu(x), (18) 

which is the Fierz- symmetric counterpart of O. We therefore extend the right hand side 
of ffTUj) to include four operators in this case: the physical operator O, the evanescent 
operators Oe and Oe> and the Fierz-evanescent operator Oe = O — O. 

The IR subtractions turn out to be particularly complicated in this case, due to 
the fact that the evanescent operator Oe already appears in the tree level calculation. 
As a consequence the naive split-up into non-factorizable diagrams, which contribute to 
the hard-scattering kernels, and factorizable diagrams, which give form factor and wave 
function corrections, is spoiled. In NLO we find that equations (30) and (31) of [5] should 
be replaced b}|§ 

+ zSHQs)® = f {0) ■ ® $ (0) + Ag, 

+ 4 1} (4} (0) = ■ T/ 0) ® m + F<°> • ff > ® - A% (19) 

where contains the (non- vanishing) 1-loop counterterms of the form factor and wave 
function corrections for the Fierz-evanescent operator Of- In other words, the split-up 
in the above example of the colour- allowed tree amplitude followed from the fact that 
the corresponding counterterms vanish for the physical operator O (i.e. A^ = 0). 

From the first equation in (fl9~l) we see that we can neglect the factorizable 1-loop 
diagrams in the computation of the NLO kernels. In order to account for the counterterm 

8 We introduce the " hat" notation to distinguish these quantities from those of the preceding section. 
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contribution A^\, we compute the UV-divergences of the 1-loop diagrams from Figure 5 
and 6 of [5] with an insertion of the Fierz-evanescent operator Op. We find that the 
counterterms (ct) are given by 



F let 



F let 



) = 2C F j\w V E {u,w) j#<°> $ (0) H + -^#1° $g } (it;) j , 



(20) 



with Ve(/u, to) from equation (45) of [5J. Convoluting these expressions with the LO 
kernels, T^p = Cp/N c and T^f = V-^c> yields the additional counterterm contributions 



AW _ r 

Lip l — L^plA F 2 



AT, 



(21) 



With this prescription the NLO kernels turn out to be free of IR-singularities. Evaluating 
the first equation of ffl~9l to 0(e 2 ) gives (in terms of T± from ffl4l ). 



T 1 (1) (w)+C F 



2N C 



2N C 

2N r . 



Cp l\L_ 

2N 2 \m\ 



2 \ e 



h(u)+2L )e+ (i 2 (u) -L 2 )e 2 + 0(e 



C h 



2 \ e 



8N* V m 
_Tf(n) 



^2 ^2L + t £)0 («)+ t B ,l(u)-i/ £ + (9(£ 



C F 
2N? 



2 \ e 



2 + fi(u)£ + £>(e 2 



:22i 



We next compute form factor and wave function corrections to O, Oe and Op. Proceeding 
along the lines outlined in Section 4.3 of [5], we obtain^ 



1 amp,E ^ E 



24 



2 

mi 



Tie 



24e(l + e) 



- 1 amp,F Y F 



-c F 
+ c F 

-Cp 

c F 



3 

- + 

e 

1 

4e ~ 



e^/i 2 V 

2 

mi 



Tie 



1 - 6e + 16e 2 



$(o) 
14e 3 



£(1 -2e)(l -£) 2 



# (0) d (0) 



4 1 



p(0) A(o) 



e^/i 2 \ £ 



1 - 



e^/z 2 V 

2 

mi 



Tie 



l-e) 

1 - 3e + 6e 2 - 6e 3 



1 - e + 2s< 



e(l - 2e)(l - e) 2 e(l - 2e) 



$(o) 



3 The corrections to the physical operator can be found in [5]. 
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2 

mi 



Tie 



1 - 3e + 6e 2 - Qe 3 



e(l-2e)(l-e) 2 



(23) 



and for the wave function corrections 



1 E ^amp,E 



48C F 
2e 



2C f 



V 



F v amp,F 



2C F 



<7f 
2e 



Vk 



^(0)^(0) 



2C F 









(24) 



where ® represents a convolution and is the Efremov-Radyushkin-Brodsky-Lepage 
(ERBL) kernel [18] (given explicitly in equation (43) of [5]). We finally compute the 
convolutions of the NLO kernels with the wave function corrections. For the physical 
operator O we get 



j?(o) fO) AO) = Lp(o) $(1) 

1 amp r> a t 2 am 



2A, 



amp 



Lp(o) t (1) + .5L 

v- 1 amp 1 2]y2 



A', 



* 3 (u) + 0(e) I F (0) $W 



(25) 



while the evanescent operators give again //-dependent corrections to the physical kernels 

1 



p(o) ^,(1) 4 (1) 
E ± l,E v amp,£ 



p(o) fO) AO) 

F ± 1,F v amp,F 



2N ( 
1 

'2JV„ 



p(0) A(i) 



p(o) ^.(l) 

-.r c J 2 F Vjnip^ 



c% 

2N* 



QL + t E)2 {u) + 0(e) 
6L + t F>0 {u) + O(e 



p(o) $(o) 



^(0) $(o) 



(26) 



According to ([12]) we now have assembled all pieces to perform the IR subtractions in 
NNLO. However, as we have seen above in the calculation of the NLO kernels, the naive 
split-up into factorizable and non-factorizable contributions is spoiled for the colour- 
suppressed amplitude. In analogy to (1191) we therefore have to account for an additional 
contribution A Fi on the right hand side of (fl2l) . which represents the 2-loop counterterms 
of the form factor and wave function corrections for the Fierz-evanescent operator Of- 

The calculation of this counterterm contribution requires a rather complicated 2-loop 
calculation on its own. We refer to Appendix [B] for the details of this calculation and 
quote the contribution to the physical kernel only, 



A 



(2) 
F,l 



C F A 



(2) 



F^F,2 



2C% 



+ 



7T 

28 + — 
12 



+ 



7L 



I ■ L)C F + Y° A ~l n f 



2 



C A n f + OU 

36 18 } v 



£(o)$(o). (27) 



3 Vertex corrections in NNLO 



As we have seen in the last section, the NNLO calculation of the hard-scattering kernels 
requires a rather complex subtraction procedure of UV- and IR-divergences. The fact 
that the kernels turn out to be free of any singularities represents both a non-trivial 
confirmation of the factorization framework and a stringent cross-check of our calculation. 

3.1 Hard-scattering kernels 

In terms of the Wilson coefficients Cj of the physical operators Qi from the operator basis 
(J3J), the topological tree amplitudes take to NNLO the form 



ai(MiM 2 ) 
a 2 {M x M 2 ) 



Air 



4vr 2N r 



Ci vl 2) + c 2 v 2 [ 



(2) 



0(al 



+ . 



a 



C 2 _^ a s C F 



N r 4vr 2N r . 



+ 



47T 



2^-|)^ + (|c 1 + |)^ + 2 c A c 2 y« 

+ 0(al)\ + ... 



+ [SC F -^C A -^n f )(' [ (' l 



(2f 



where the ellipsis refer to the terms from spectator scattering which we disregard in 
the following. In this notation the a s corrections have been expressed in terms of the 
convolution 



J du ^ — 6L + g 2 (u) + iirg^ufj 



(29) 



where L = In /z 2 jm\ and (recall that u = 1 — u 
gi (u) = — 3 — 2 In m + 2 In m, 
3(1 -2m) 



02 W 



-22 + 



lnii + 



u 



2Li 2 (u) — In 2 u — In u — (u — > u) 



u 



(30) 



If we transform these expressions into the Fierz-symmetric operator basis that has been 
used in many previous QCD factorization analyses, we reproduce the NLO result from pQ. 
In NNLO we find the convolutions 



Vl 2) = f^ du { ( 36C V - 29C A + 2t2/)l 2 

f /29^ 2 \ . , 91 10 



+ m 



/29 2 



L 



10 



+ C F h 7 (u) + C A h 8 (u) + (n f - 2)h 9 (u] 0) + h 9 (u; z) + h 9 (u; 1) 



+ in 



C F h 2 (u) + CAh 3 (u) + (n f - 2)hi(u; 0) + h A (u\ z) + /i 4 (w; 1) 



t>M 2 {u) 



V 2 {2) = ! du |l8L 2 + (^21 -6g 2 (u) - Ging^uj^L + h 5 {u) +inho(u)\<f> Ma (u), (31) 





where n/ = 5 represents the number of active quark flavours and z = m c /mb. The 
explicit expressions for the NNLO kernels h _4, which specify the imaginary parts of the 
topological tree amplitudes, can be found in [5]. As a new result we obtained the real 
parts of the topological tree amplitudes to NNLO, which have been given in terms of a 
new set of kernels h^g that are listed in Appendix O 

Partial structures of our NNLO result can be cross-checked. First, we verified that 
the scale dependence between the Wilson coefficients, the coupling constant, the hard- 
scattering kernels and the light-cone distribution amplitude cancels in the tree amplitudes 
a i (M 1 M 2 ) to 0(a 2 s ) as it should^. Second, we compared the terms proportional to rij 
with the analysis of the large /3o-limit in [TH] and found agreement. Finally, we reproduced 
the imaginary part of the colour-suppressed amplitude from our earlier analysis in [5], 
which was derived on the basis of Fierz-symmetry arguments. 



3.2 Convolutions in Gegenbauer expansion 



We expand the light-cone distribution amplitude of the emitted meson M 2 into the eigen- 
functions of the 1-loop evolution kernel, 



>M 2 W 



Quu 



n=l 



(32) 



where a^f 2 and Ci 3 ^ are the Gegenbauer moments and polynomials, respectively. It is 
convenient to truncate this expansion at n — 2, which allows us to perform the convolution 
integrals in our final expression (|3T|) explicitly. The convolution with the NLO kernel 
results i 



f 1 , ( s, , s 45 11 M2 21 
J du g 2 {u) (pM 2 (u) = ~y + y a i ^ a 



M 2 
2 5 



(33) 



whereas the convolutions with the NNLO kernels become 



du h 5 (u) 4>m 2 (u) 



5347 14833 



60 



<3 + 3744C 5 



12487 
12 



936Cs + 72 In 2 W 2 + 



239vr 
90 



10 We emphasize that this cancellation would have been incomplete, if //-dependent contributions from 
the mixing of evanescent operators as e.g. in (|16p or (|26[) had been missed. 
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We refer to [5] for the convolutions with the kernels g\ and /io-4- 
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,'4568 77157, , rt „ rtn> 



f 32369221 2236872 _ 
+ 1^2600 3^ C3 + 74304 Cs 



+ 



18576C 3 - 20641n2)7r 2 + 



/' 

Jo 



du he(u) M » = 348 - ^ af 2 + — af 2 , 



204218 u _\ , 797tt 4 

9 

154 „„ 329 



10 



7 M 2 
*2 5 



1 , , , , , / \ 12809 26606 _ 
du h 7 {u) (pM 2 {u) = — ^ — C3 + 6564 C 5 



+ (Hli- 176 4 C3 -48ln2)^ + 13414 



45 



, 50387 132294 im _ 

+ { + Cs - 32472 Cs 



/66425 M , \ , 176tt 4 1 
( — - 8856 C 3 - 1296 In 2 J tt 2 — \ 



[75807647 3960924 , nftrtftil/ . 

+ < Cs + 129204C 5 

\ 12600 35 S3 S5 



/2074841 n , nn , A n „„ rtl n \ 2 727tt 4 1 
+ [ — — 34884 Cs - 8672 In 2 J tt 2 + — ^— ^ a 



M 2 
2 5 



JO 



, , , x , /x 74611 618 A 

du h 8 {u) (pM 2 {u) = — — + — C3 - 186 C 5 

180 5 



_(5i5_ 10 8C 3 -361n 2 )^- 169 ^ 



120 



[355693 10818 

+ {^ + ~ C3 " 2556(5 



1081 t \ o 1517T 
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270C 3 + 6841n2)7r 2 + 



,M 2 



148920211 128283 , 



66545 _\ 2 1403tt 4 , ;/ 



18 1458 C3-41201n2). 



(34) 
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We finally perform the convolution with the hard-scattering kernel hg(u; Zf), which stems 
from the diagrams with a closed fermion loop. As this contribution depends on the mass 
rrif = Zfrrib of the internal quark, we parameterize the convolution as 
-i 



du h 9 (u; z f ) 4>m 2 (u) = H 9<0 (zf) + #9,1(2:/) af 2 + H 9j2 (zf) 



M 2 



For massless quarks we may perform the convolution integral analytically, 



du h g (u;0) 4>m 2 {u) 



493 2vr 2 



40 



+ 2tt 2 ) af 2 



8059 
~60Y 



(35) 



(36) 



whereas we obtain numerical results for massive internal quarks. In Table[T]we summarize 
the contributions from closed fermion loops for massless quarks (z q = 0), for a 6-quark 
(z b = 1) and for a charm quark [z c G [0.25,0.35]). 

We illustrate the relative importance of the individual contributions setting fi = nib 
and z c = rric/mt, = 0.3, which yields (with Cf = 4/3, Ca — 3, n/ = 5) 

yW = (-22.500 - 9.425 i) + (5.500 - 9.425 i) af 2 + (-1.050) af 2 , 



V} 2) = (-178.38 - 349.44 i) + (660.59 - 119.36 i) af 2 + (-85.40 - 62.63 z) a Ma 



2 ) 



y 2 (2) = (322.19 + 320.94 i) 



;-212.97 + 154.41 i) af 2 



(3.81 - 34.06 i) a 



M 2 
2 ■ 



(37) 



We find relatively large coefficients for the NNLO terms and expect only a minor impact 
of the higher Gegenbauer moments in the symmetric case with af 2 = 0. 

We conclude with a remark concerning the large /?o - h m it that has been considered 
in [12]. In this approximation we get 



(2)1 



I 00 



-239.31 - 264.94i) + (380.33 - 252.90i)af 2 + (-40.96 - 21.68z>f 2 , (38 



,m 2 

l 2 > 



whereas the contribution from V 2 is completely missed. As a consequence the NNLO 
contribution to «i is substantially underestimated in this approximation, whereas the 
one to a?2 deviates from the full NNLO result between ~ 15% for the imaginary part 
and ~ 40% for the real part. This illustrates the importance of performing exact 2-loop 
calculations. 



z f 





0.25 


0.275 


0.3 


0.325 


0.35 


1 




20.81 


17.12 


16.43 


15.72 


14.99 


14.26 


-3.62 




-33.07 


-13.28 


-12.37 


-11.54 


-10.77 


-10.07 


-0.68 


Hg,2 


3.56 


2.08 


1.94 


1.81 


1.68 


1.57 


0.01 



Table 1: Fermionic contribution in the notation of $35\) . The first column 
refers to massless quarks, the last column to the b-quark and the other columns 
to the charm quark for different physical values of z c = m c /mb. 
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4 Numerical analysis 



We conclude with a brief analysis of the numerical impact of the considered NNLO cor- 
rections. As a phenomenological analysis of hadronic B decays is beyond the scope of the 
present paper, we focus on the perturbative structure of the topological tree amplitudes 
and discuss their remnant uncertainties. In particular, we now combine our results with 
the NNLO corrections from 1-loop spectator scattering that have been worked out in [3]. 

4.1 Implementation of spectator scattering 

In contrast to the vertex corrections considered in this work, the spectator scattering 
term is sensitive to two perturbative scales: the hard scale fih ~ Tnj, and a dynamically 
generated intermediate (hard-collinear) scale Hh c ~ (^QCD m b) 1 ^ 2 ■ The hard scattering 
kernels from spectator scattering therefore factorize further into coefficient functions H* 1 , 
encoding the hard effects, and a universal hard-collinear jet-function Jn. Renormaliza- 
tion group techniques can be used to resum parametrically large logarithms of the form 
lnmfe/AQCD in terms of an evolution kernel tl\\. Following the first paper of [3], we 
implement the spectator scattering contribution to the topological tree amplitudes ag[jf| 

-» Ci{ph) H^ijlh) ®U\\(Hh,H hc ) <g> J\\(lMhc) ® [fB<j>B](phc) ® <pNh{^hc) ® 4>M 2 {^h)- (39) 

Since the spectator scattering starts at 0(a s ), the resummation is required here in 
the next-to-leading-logarithmic (NLL) approximation. Unfortunately, a complete NLL 
resummation is not possible since the evolution kernel Un is known in the leading- 
logarithmic (LL) approximation only [2"U] . 

We therefore proceed along the lines of our earlier analysis where we worked in 
the LL approximation which is consistent for the imaginary parts that are of 0(a 2 s ). 
According to this, we implement the LL evolution of the HQET decay constant and the 
Gegenbauer moments to evolve the hadronic parameters from their input scales to the 
ones required in (I3"9l . The B meson distribution amplitude is modeled according to [21J, 
which implies Ab(IGcV) = (0.48 ± 0.12)GeV and, for the first two logarithmic moments, 
eri(lGeV) = 1.6 ± 0.2 and <72(lGeV) = 3.3 ± 0.8. The 1-loop matching corrections to 
the hard functions H- 1 [3] and the jet function Jn [20j [22] are implemented neglecting 
crossed terms of O(a^). We finally adopt the BBNS model from pQ to estimate the size 
of power corrections to the factorization formula. 

In the spectator scattering term we compute the Wilson coefficients from the effective 
weak Hamiltonian in the NLL approximation with 2-loop running coupling constant. 
Quantities referring to the hard scale are evaluated in a theory with n/ = 5 flavours and 
those referring to the hard-collinear scale with rif = 4. 

12 One should keep in mind that the Wilson coefficients in the spectator scattering term refer to a 
different operator basis than the one used in the current work (namely the Fierz-symmetric traditional 
basis that we denoted by a tilde in [5]). 
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4.2 Tree amplitudes in NNLO 



We finally evaluate the topological tree amplitudes for the B — > 7777 channels using the 
input parameters from our earlier analysis [5] and computing the Wilson coefficients 
in the vertex corrections in the next-to-next-to-leading-logarithmic (NNLL) approxima- 
tion [SI E3] with 3-loop running coupling constant [24] and = 205 MeV. Under these 
specifications the NNLO prediction of the topological tree amplitudes become J^l 

01(7777) = 1.008| y(0) + [0.022 + 0.009i] yW + [0.024 + 0.026?:] y(2) 
- 0.012| s(1) - [0.014 + 0.011i] 5(2) - 0.007| p 

= 1.019±gSI + (0.025±gSi)i, 
a 2 (7r7r) = 0.224| y{0) - [0.174 + 0.075z] - [0.030 + 0.048.] y(2) 
+ 0.075| s(1) + [0.032 + 0.019z] s(2) + 0.045| 



\p 



= 0.173™ -(0.103^)*. (40) 

Here we disentangled the contributions of the various terms in the factorization formula, 
namely the tree level result V^ ' ("naive factorization"), NLO (1-loop) vertex corrections 

NNLO (2-loop) vertex corrections V^- 2 \ NLO (tree level) spectator scattering 
NNLO (1-loop) spectator scattering and the modelled power corrections P. 

The new contributions from this work consist in the real parts of the terms denoted by 
y( 2 ) _ p or colour-allowed amplitude 0:1(7777), this correction is slightly larger than the 
a s terms due to an numerical enhancement from the Wilson coefficients in the effective 
HamiltomarH On the other haud, the colour-suppressed amplitude o 2 (™) receives a 



13 The numbers for the imaginary parts differ slightly from those of 5], since we now evaluate the 
Wilson coefficients throughout in the NNLL approximation. 

14 We remark that a similar enhancement is unlikely to exist at even higher order of the perturbative 
expansion, since the NNLO expressions already reveal the full complexity. 







Hhc 


f B 




Xb 




X H 


Re(ai) 


+0.008 
-O.Off 


+0.006 
-0.007 


+0.003 

-0.003 


+0.006 
-0.008 


+0.006 
-0.009 


+0.007 
-0.008 


+0.007 
-0.007 


Im(a!i) 


+0.017 
-0.011 


+0.002 
-0.003 


+0.001 
-0.001 


+0.002 
-0.003 


+0.002 
-0.003 


+0.004 
-0.004 


+0.007 
-0.007 


Re(a 2 ) 


+0.016 
-0.008 


+0.026 
-0.023 


+0.014 
-0.014 


+0.038 
-0.025 


+0.039 
-0.026 


+0.038 
-0.033 


+0.045 
-0.045 


Im(a2) 


+0.019 
-0.028 


+0.005 
-0.004 


+0.002 
-0.002 


+0.005 
-0.003 


+0.005 
-0.003 


+0.007 
-0.006 


+0.045 
-0.045 



Table 2: Dominant uncertainties of our final predictions for the colour- allowed 
tree amplitudes 01(7777) and the colour-suppressed tree amplitude 02(7777) from 
scale variations, hadronic input parameters and modelled power corrections. 
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moderate correction. In particular, we do not find an enhancement of the phenomeno- 
logically interesting ratio |«2/«i| from the perturbative calculation. 

In Table [2] we list the uncertainties of our NNLO predictions stemming from scale 
variations, hadronic input parameters and the modelled power corrections. The values 
of the first two columns follow from varying the perturbative scales independently in 
the ranges [L^ = 4.8l 2 !4 GeV and ^hc = l-5lo!s GeV. As the dependence on the hard 
scale tends to cancel between vertex corrections and spectator scattering, we vary both 
contributions independently and take the larger interval (from the vertex corrections) as 
our estimate for higher order perturbative corrections. The scale dependence of the vertex 
corrections is also illustrated in Figure [U where we read off that it gets substantially 
reduced for the real parts at NNLO, whereas the reduction is less pronounced for the 
imaginary parts. 

For our final error estimate in (HOl) we added the individual uncertainties from Ta- 
ble |2] in quadrature. Whereas the colour-allowed amplitude ai(7T7r) can be computed 
precisely in the factorization framework, the situation is less fortunate for the colour- 
suppressed amplitude a 2 (7r7r). Due to large cancellations between the vertex corrections, 
the colour-suppressed amplitude becomes particularly sensitive to the spectator scatter- 
ing contribution and is therefore subject to rather large uncertainties related mainly to 
our restricted knowledge of the hadronic input parameters. 




Im(ai) 




Re((W 



v 



/<■/, 



Im(a 2 ) 




Figure 2: Dependence of the tree amplitudes ai(inr) as a function of the hard scale fj,h (vertex 
corrections only). The dotted (black) lines refer to LO, the dashed (orange/light gray) lines to 
NLO and the solid (blue/dark gray) lines to NNLO. 
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5 Conclusion 



We computed the real parts of the 2-loop vertex corrections for charmless hadronic B 
meson decays, completing the NNLO calculation of the topological tree amplitudes in 
the QCD factorization framework. We in particular showed how to compute the colour- 
suppressed tree amplitude without making use of Fierz-symmetry arguments and found 
that the hard-scattering kernels are free of IR-singularities and the resulting convolutions 
with the light-cone distribution amplitude of the emitted light meson are finite, which 
demonstrates factorization at the 2-loop order. 

The numerical impact of the considered corrections was found to be moderate, al- 
though they can be of similar size as the NLO corrections. The scale dependence of 
the real parts of the topological tree amplitudes is significantly reduced at NNLO, which 
allows for a precise determination of the colour-allowed amplitude a±. In contrast to this, 
it remains difficult to compute the colour-suppressed amplitude a 2 in the factorization 
framework, since it is subject to substantial uncertainties from hadronic input parame- 
ters and potential l/ra& corrections. In particular, we do not find an enhancement of the 
phenomenologically important ratio |a 2 /«i| from the perturbative calculation. 
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A Auxiliary coefficient functions 

In the calculation of the colour-allowed tree amplitude, the NLO kernels have been given 
in (TH|) in terms of the coefficient functions 



lnw lnu\ 7T 2 

22, 



t (u) = 4Li 2 (w) - ln 2 u + 21nwlnu + ln 2 w + (2 - 3m) f 

V u u J 3 

t x (u) = -2Li 3 (u) - 2Si )2 (w) - 2]nuLi 2 (u) + In 3 w - 21n 2 wlnu + lnwln 2 w - In 3 w 

2-3u 2 T .,, 2-3u/, , , , _\ 6-11m + 2m7t 2 , 
H = — Li 2 (u) — In u — In u In u ) H In u 

uu u V / u 

4 -3m, 2 18 - 33u + 5mr 2 , (7 - 6u)n 2 

+ — ln 2 n \nu+ K - + 2( 3 - 52, 

2u 6u bu 

t 2 (u) = 10Li 4 (w) - 8S 2i2 (w) + 10Si, 3 (u) -81ntZLi 3 (u) + 101ntZSi, 2 (w) - ^ln 4 w 

4 17 

+ 5 In 2 u Li 2 (u) H — In 3 u In u — In 2 u In 2 u H — In u In 3 u H In 4 u 

v ; 3 3 12 

2-6m + 6m 2 t . , s 4-6m + 3m 2 / m , N , , _ T . , x \ 8 - 3m , 3 _ 



H z Li 3 (w) = (Si 2 (u) +\nuLi 2 (u) ) In u 

uu uu V / bit 



17 



2 - 3m / , o , o , , , 2 \ 60 1 - 2m + 17mtt 2 , 2 

H 41n 3 M - 61n 2 MlnM + 3hiMln 2 M \n 2 u 

6u V / 12m 

3(6 - 4m - 7m 2 ) + mmtt 2 , , 24 - 54u + 5utt 2 , (29 - 24m) 

+ — Li 2 (m + — InM In u + — - 

Suu 6m 6m 

6(12 - 13m) + 7mtt 2 2 _ 24(7 - 13m) + (10 - 15m)vt 2 23tt 4 
+ 12m 11 U+ 12M nU ~ 180 

24m(7 - 13m) + (2 + 23u - 27u 2 W 2 + 24mm£ 3 , _ 10- 11m > „ _ 

InM H = (3 - 112, 



7T 



12mm m 
1 — 2m /InM lnM\ 16 



1 — 2m t . , . 1 — 3m l2 u , 1- 2 — 3m , 2 - 
= — Li 2 (mJ H — In u -\ — -mMlnM m u 

2mm 4m 2m 4m 

4(1 — 2m) /InM lnM\ (6 — 5m)7T 2 



, 12 (41) 
V u u J 12m v ; 

and the convolutions of the NLO kernels with the wave function corrections, cf. (fT5"j) and 
(flEj) . involve 

* 3 (u) = 4Li 3 (M) + 4S 1i2 (m) - 4 InM Li 2 (u) + ^ In 3 m - 2 In 2 m InM - \ In 3 M L ' 2 ' " ' 



3 3 MM 

1-3m/ , 2 n _, , _ _ n 2 -\ 3 , (4-3m)tt 2 15 

Mm m + 2u In m in u — u In m in u H 4f 3 , 

2um V J 2m 6u 2 S3 ' 



t 4 (u) = 12Li 4 (M) - 20S 2j2 (m) + 12S 1j3 (m) - 8 ^ In m + In M^ Li 3 (m) + 121nuS 1)3 (u) 

+ 41nuSi j2 (u) + ^41n 2 M + 41nMlnu + 21n 2 ujLi 2 (M) - ^ln 4 u+ ^ln 3 Mlnu 



1,2,2- 1, ,3- 3 l4 _ 4-llu + 3u 2 T . . . 5 -12m, 3 
In Mm m Inuln mH — in u Li 3 m H In u 

2 3 4 mm 6m 

1 + u-3m 2 , , 2 - 10m + 6m 2 , t . , , 1-5m + 3m 2 , _ t . /n 
H bi 2 (m) H = lnMLi 2 (M) lnMLi 2 (M) 

UM MM UU 

2 - 10m + 9u 2 , 2 , _ 1 - 2m , , 2 _ 5 - 6m , 3 _ 
H = In m In m — In u In u In u 

2mm 2mm 6m 

18 - 24m + 15u 2 - lOumr 2 , , 16 - 27u + 4mr 2 , 2 

= Li 2 (m) — In m 

3mm Au 

6 - 36m + 27m 2 - Auuix 2 , , _ 3(14 - 17 u) + 8wr 2 , , _ 

= In u In it H In u 

2mm 12m 

8 - 15m - 4tt 2 - 48m( 3 3(2 - 3m)C 3 23vr 4 (23 - 17m)tt 2 

+ 4M nU+ M W + 12M 

81 - 126u + 45u 2 - (14- 22u + 6u 2 )tt 2 - 192u< 3 _ 137 



12mm 



InM 
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6(1 — 2u) 6 7r 2 

t E2 (u) = — = — -Li 2 (w) \nu\nu - 6\nu - 6\nu - — + 50. (42) 

uu u u 

In the calculation of the colour-suppressed tree amplitude, the NLO kernels in f[2"2"j) contain 
the coefficient functions 

it 

ti(u) = — In u — In u + 8 + iir, 
u 



r 2 

20 



" / \ U I ' t ■ i \ i 2 i i - a i \ 1, 2- - (3 — 2u)7T 

t2{u) = — Li 2 (m) — In u + In u In it + 4 In u I H — In w — 4 In w = h 

u V / 2 6m 

+ 27r^4 — lnwj , 

w 

^£ o( M ) = — ln« — In w + 6 + Z7T, 
u 

t E ,i{u) = --(Li 2 (u) + ln 2 u-3\nuj + ^ In 2 it - 3 lnw + 14 - y + m (3 - \nuj (43) 

and the convolutions with the NLO kernels, cf. fT25|) and ff26l) . give rise to 
Ti" 2 _ „, ,o_ hi w 



i 3 (u) = 5 + — In u + 2 In u In w — In w 

3 m u 

t B2 (w) = -(6Li 2 (w) + 3wlnw-7r 2N ) - 3 -^-^^ In w + 27 + 3ivr, 
it V J u 

tFoM — — (^2(1 + 2-u)Li2(w) — wln 2 -u — 2-uln-uln-u + 3-uln-u — (2 + u) — 
it V 3 / 

+ *] Jl *u- 3(1 + ^ lnn + 29 + ^3- — lni/+ — In A (44) 

11 U V U 1i / 

B Calculation of 2-loop counterterms A^ 

We present the calculation of the 2-loop counterterms A^, that are required in the 
NNLO calculation of the colour-suppressed tree amplitude as described in Section 12.41 
The counterterms receive three contributions 

Ag = f$ ® ji^ |f + PP + PP } ^ (45) 

where <S> represents a convolution and "ct" refers to the counterterm contributions of the 
form factor and the wave function corrections. Notice that the wave function corrections 
actually correspond to local corrections to the decay constant, as a consequence of the 
fact that the tree level kernels T-°p are constant (cf. also (T201) and (T2TT) ) . 

We first consider the mixed term Fp &p , which involves the calculation of the dia- 
grams from Figure [31 The first diagram vanishes due to a scaleless loop integral and the 
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Figure 3: Diagrams that contribute to the mixed contribution F F & F . The 
symbol in the lower (upper) line refers to an insertion of the 1-loop counter- 
term from the form factor ( wave function ) correction of the operator O p . 



second diagram yields 

1 _ 6e + 16e 2 - Ue 3 ~ (0) ^ (0) 1 ^(o) £(o) \ U(]] 

+ 4^(l-2e)(l-^ Fe ® e + 16.(1 - ef ^' ® E ' )' (46) 

We are left with the 2-loop counterterm from the last diagram of Figure El which requires 
the calculation of the UV-divergences of the 2-loop diagram from Figure HI For this it is 
convenient to apply the method proposed in [25J (sometimes called IR-rearrangement) , 
which allows to set all masses and external momenta to zero. The calculation then reduces 
to the evaluation of 2-loop tadpole integrals, which depend on a single mass scale (an 
artificial scale that has been introduced to separate UV- and IR-divergences) . Computing 
the 1-loop counterterms with the same prescription and accounting for the wave-function 
renormalization, we get 

= + 5) ^ (0) $ (0) - ( - l) FP *<? + J^ F e> + -/P (47) 

Next we compute the form factor correction F F 2 ^ (the corresponding diagrams are 
shown in Figure [5^i). The first diagram gives again the contribution S± from fj46|) . On 
the other hand the computation of the 2-loop counterterm from the second diagram 
of Figure Eh is rather involved. It requires the calculation of the UV-divergences of a 
couple of 2-loop diagrams (shown e.g. in Figure 1 of [15J) and the corresponding 1-loop 
counterterms. Proceeding as before with the method of IR-rearrangement and accounting 
for the 2-loop wave-function renormalization in the MS-scheme [26J, 



zfl = zg = c F \ ( \c F + c A ) I + (\c r - H c , + 1„,) (4£ 




Figure 4: The UV-divergences of this 2-loop diagram contribute to 62- 
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Figure 5: Diagrams that contribute to 6^ (a) and F^ 1 <E>^ ; (b) 



yields the 2-loop form factor counterterm for the Fierz-evanescent operator O 
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32e : 



+ 



64 32 e 
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6 : C , + _ n/j _^ 4 <»» 



(49) 



We finally account for the wave function correction F^ -* from Figure Eb- The first 
diagram again vanishes due to a scaleless integral and the second diagram yields, in a 
convolution with a constant kernel, again the contribution 83 from (T49|) . 

To summarize, in terms of the individual contributions 8i from (|46p . (T4"7]) and (|49p . 
the 2-loop counterterms required in the calculation of the colour-suppressed amplitude 
become 



Agi = C F A% = ^ 28 1 + 8 2 + 28 3 



(50) 



C NNLO hard-scattering kernels 



Our final expressions for the real parts of the NNLO vertex corrections from ( 13~TT) involve 
the following set of hard-scattering kernels, 



4(3 - 3m + 8m 2 - 2m 3 ) T . . 8(12 - 35m + 36m 2 - 14m 3 + 4m 4 ) 

Li 4 (m) ^ S 2 , 2 (m) 



M 3 M 2 
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+ 
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5— Li 2 (m) 2 + 41n 2 MLi 2 (M) + - ln 4 M ln 3 MlnM 

2u 6 u A 3 3 

122 - 91m 2 + 20m 4 



3mm 3 



hiM hiMLi 2 (M) — 



51 + 16m 2 + 268m 4 + 4m 6 
12m 3 m 3 



In Mln m 



15 We performed this calculation for arbitrary bilinear quark currents, which allows us to perform 
several cross-checks. We in particular verified that the anomalous dimension of the vector current 
vanishes at the 2-loop level and reproduced the one of the scalar and the tensor current from [27] , 
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M 3 M 3 



8 hiM L1 3 (m) + In Mln u — 8(^3 hi m 
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2(3 - 11m + 15m 2 - 8m 3 + 4m 4 ) / , _ 7tt 2 \ /n 

' InMlnM Li 2 (m) 
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3 - 102m + 293m 2 - 278m 3 + 36m 4 + 74m 5 - 43m 6 + 14m 7 T . . 
+ — Li 3 M 

MM 
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6m 2m 3 m 
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+ ~ 3 ( Li 3 (-m) -lnMLi 2 (-M) ln(l + M)J 



22 



+ 
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21 - 62m - 29m 2 + 182m 3 - 91m 4 



Li 2 («) + 



42 + u - 142m 2 + 91m 3 



1ii 2 m 



2m 2 m 2 ' 4mm 2 
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12m 3 m 3 
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IiimIiim 



12m 2 m 
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InM 
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90m 3 m 3 

(m - m)(3 - 110m + 96m 2 + 28m 3 - 14m 4 ) ( 3 



+ 



24m 2 m 2 



h 6 (u) 



2m 3 m 3 

327 3(1 - 2m) 2 3(1 -2m 2 ) 



— (m <-> u) 



+ 



2 2m 
(1-2m 2 )tt 2 



2mm 



In u In u — 



3(13 - 24m) 
2iZ 



InM 



4mm 



+ (m <-> u) 



+ 



8Li 3 (u) — 8 In u Li 2 (u) + - In 3 u — 4 In 2 u In u = Li 2 (u) 

3 MM 



25 - 24m, 2 13, , 9 , lbr 2 . 
H = In m H — — In m In m m u = (u <-> m) 

2m m 2m 6m 



^7(m) 



(1+M)( 3^f + 3 " 2) (l2«i(«) + vr 2 ln(l + „)) - 48 InM U 3 (u) 



+ ^(24H 2 (M)-27T 2 Li 2 (-M) ) + 



4(6 - 9m + 16m 2 - 4m 3 ) 



Li 4 (M) + 81n mL1 2 (m) 



4(52 - 152m + 156m 2 - 61m 3 + 18m 4 -m 5 ) m . . 4, 4 8,, , 

^S 2j2 (m) + -ln 4 M - -ln 3 MlnM 

o o 



M 3 M 2 



4(52 - 204m + 308m 2 - 209m 3 + 107m 4 - 3m 5 + M e 

M 3 M 3 



lnMLi 3 (M) — ^ 3 InM 



13 - 54m + 88m z - 84m 3 + 82m 4 - 48m 5 + 16m 6 t . 2 3 - 18m + 6m 2 - 4m 3 3 



M 3 M 3 
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M 3 M 3 



-Li 2 (M) 2 + 



InM lnMLi 2 (M) 



6m 



ln^M 
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M 3 M 



InM Li 2 (u) 
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2m 3 m 3 



In m In u 
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f — ^Li 3 (M) 
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6 - 116m + 149m 2 - 50m 3 - 6m 4 + 6m 5 -4m 6 , 2 , 
— In m In m 

2m ,4 m 

2(3 - 3m - 7m 2 - 3m 3 ) / . . , . . In 2 u + tt 2 , 
+ — § Li 3 (-m) - lnMLi 2 (-M) ln(l + M) ) + (m - u) 

XL \ Zi 



/92 - 117m + 109m 2 + 16m 3 - 8m 4 (1 - 2m + 2m 2 )(83 - 14m + 14m 2 )tt 2 \ t . . 
X { W + 6m^ J Ll2(M) 

/ 46 - 73m + 26m 2 + 13m 3 -4m 4 , 2 Y , / 92 - 289m + 421m 2 - 264m 3 + 132m 4 
+( 7T—t 47r 2 )ln 2 M+( - 



\ 2mm 2 / \ 4mm 2 

(139 - 554m + 856m 2 - 558m 3 + 164m 4 + 138m 5 - 46m 6 )tt 2 \ . 

H 9 _o ) InMlnM 

3m^m j 



InM 



+ 



+ 



183 - 308m (96 - 160m + 55m 2 + 15m 3 + 3m 4 - 2m 5 )tt 2 
2m 3m 2 m 

(191 - 784m + 1262m 2 - 1121m 3 + 973m 4 - 495m 5 + 165m 6 )tt 4 

180m 3 m 3 

(580 - 1763m + 2239m 2 - 952m 3 + 476m 4 )tt 2 
48m 2 m 2 

3(1 - 37m + 119m 2 - 138m 3 + 4m 4 + 78m 5 -26m 6 )C 3 1659 , 
^ 1 a h (m ^ m) 

MM 4 

(1 +m)(1 + 4m - 7m : 



3mm 2 



12Wi(«) + tt 2 ln(l + u) - 36 InM U 3 (u) + 14 In 2 u Li 2 (u) 



2(1 -4m + 3m 2 -m 3 ) / ini . , n 2t . . A 4(2 - 23m + 18m 2 - 6m 3 ) T . . , 
+ — 3 ( 2m 2 {u) - 27r 2 Li 2 (-M) j + -i — J -Lk{u) 

4(52 - 152m + 156m 2 - 61m 3 + 18m 4 -m 5 )^ . . (u - m)(13 - 22mm) t . x2 
- — 2 -S 2)2 («) - - ^Vs -Li 2 « 

MM MM 



2(13 - 51m + 77m 2 - 59m 3 + 13m 4 - 7m 5 + m 6 ) 



M 3 M 3 



81nMLi 3 (M) + In Mln u — 8( 3 \nu 
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— In m In m Li 2 (m) — In u H — In ' 

mm V 6 / 3 



2(3 - 110m + 339m 2 - 362m 3 + 100m 4 + 59m 5 - 46m 6 + 14m 7 ) t . , 
+ - ^ -Li 3 (m) 

6 - 156m + 227m 2 - 57m 3 + 38m 4 , T . , , 9 - 32m + 14m 2 , 3 
— In m Li 2 [uj = In m 

MM 6M 



6 - 104m + 165m 2 - 54m 3 + 34m 4 - 14m 5 
2m 3 m 



In 2 m In m 



2(3 + 3m + 5m 2 + 3m 3 ) / . . . , . . ln 2 M + 7r 2 , . A 
+ — ~ 3 -f Li 3 (-u) -ln«Li 2 (-u) ln(l + M)J 
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/92 - 381m + 161m 2 + 440m 3 - 220m 4 4(1 - 3m + 3m 2 )(m + 2m 3 - m 4 )tt 2 \ 
+ I 4W 3m%5 J Ll2(M) 
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2mm 2 3 7 90m 3 m 3 

/(l +m) (246- 169m 2 ) (87- 210m 2 - 16m 4 )tt 2 ' 
\ 4m 2 m 2 3m 2 m 3 

/93 (96 - 112m - 31m 2 + 34m 3 - 7m 4 )tt 2 | , ^ 
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Au 6 u 
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/4 - 3m + 24m 2 - 42m 3 + 21m 4 (127 - 590m 2 + 652m 4 - 128m 6 )tt 2 \ , , 
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V 8m 2 m 2 48m 3 m 3 / 

/1433 - 2710m (12 - 30m + 8m 2 + 11m 3 + 3m 4 - 2m 5 )tt 2 \ , 
+ { 24m" 6m% ) lnU 
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360m 3 m 3 

(6 - 57m + 134m 2 + 26m 3 - 463m 4 + 540m 5 - 180m 6 )C 3 12641 , 

H TT^ To — + (m ^ m) 

4m 3 m 3 48 

(1+m)(1 6 ^ M ~ 7m2) (lZHi (m) + 7T 2 ln(l + ufj + 2 In u Li 3 (u) - 3 In 2 m Li 
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H — I 8mMLi 3 (M) + in mIii it-eQJnw 
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H n _ n I InM InM — I Li 2 (m) — — m u — — m umu 
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12m 3 m 
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[u - u) (3 - 10m + 10m 2 ) (33 - 21m + 35m 2 - 28m 3 + 14m 4 )tt 2 



6m 2 m 2 12m 2 m 3 
1009 (12 + 26w - 121m 2 + 80m 3 - 7m 4 )tt 2 



IiimIiim 



+ 4C 3 InM 



+ 



72m 6m 2 m 
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2880m 3 m 3 
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(51) 



The definition of the functions Hi }2 (x) can be found in Section I2~2l The diagrams with 
a closed fermion loop give for massless internal quarks 



hg{u] 0) 



12 

(1 + m)tt 2 



125 Li 2 (u) l-3u, 2 1 + u, , _ 17(1 - 2m) , 
H 1 In m H — In u In u z In u 



u 



+ (m <-> u) 



2u 



2m 



Qu 



+ 



12m 

4 T . . . 2 , 3 4 . 2 . _ 32 - 29m t ... 35 - 29m , 2 
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3 3 3 9m 18m 



3m 



In m In m 



13 + 24mV 

vm 



InM + 



7T 



18m 



m «-> u 



(52) 



and for an internal 6-quark 



h 9 (u; 1) 
1 



— Li 3 (-x b ) - Si, 2 (-x b ) - ln(l + x b )\A 2 (-x b ) - — In 



1 1.3 %b 



rr 



x b 
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6 b I + x b 
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9m 3 m 3 
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u 
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^ 4m 



Li 2 (-x b ) - - In 2 + I ln 2 (l + Xb ) + ^ ) + ^ L1 3 (m 

4 1 + x b 2 6 
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L1 2 (m) 



w 
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+ 



24mm 
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2y b (38 + 29m) f , 1 2 1,2^, x, ^ 

+ o Li 2 -x 6 - - In — — + - In (1 + x h ) + — 

9m \ 4 1 + Xb 2 6 

4(1 + 3m 2 - m 3 ) T . 128 - 504m + 389m 2 



Li 3 (w) — lnw 

ii. 6 . . 
Li 2 (u) 



3m 3 " 18m 2 m 

32 - 204m + 504m 2 - 584m 3 + 405u 4 - 150m 5 + 29m 6 



9m% 3 

+ (42 - 49u + 39m 3 )tt 2 + 4( 3 + 261 - 325m 2 
54mm 3 m 2 m 3 9mm 2 

where we introduced the shorthand notation 



(53) 



x b = ^(yb - 1), yb = y—^—- ( 54 ) 

We finally refrain from presenting the charm quark contribution, which is rather com- 
plicated and depends on two parameterizations for the 4-topology Master Integrals that 
we could not solve in a closed analytical form (cf. the discussion in [15]). We may still 
evaluate these Master Integrals numerically in Section 13.21 to perform the convolution 
with the light-cone distribution amplitude of the emitted meson M%. 



References 

[1] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Phys. Rev. Lett. 83 
(1999) 1914 [arXiv:hep-ph/9905312] ; 

M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nucl. Phys. B 591 (2000) 
313 |arXiv:hep-ph/0006124] ; 

M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nucl. Phys. B 606 (2001) 
245 |arXiv:hep-ph/010411 0|. 

[2] C. W. Bauer, S. Fleming, D. Pirjol and I. W. Stewart, Phys. Rev. D 63 (2001) 
114020 [arXiv:hep-ph/0011336] ; 

C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D 65 (2002) 054022 
|arXiv:hep-ph/0109045] ; 

M. Beneke, A. P. Chapovsky, M. Diehl and T. Feldmann, Nucl. Phys. B 643 (2002) 
431 |arXiv:hep-ph/0206152]. 



[3] M. Beneke and S. Jager, Nucl. Phys. B 751 (2006) 160 |arXiv:hep-ph70 512351|; 
N. Kivel, JHEP 0705 (2007) 019 larXiv:hep-ph/0608291| ; 
V. Pilipp, PhD thesis, LMU Miinchen, 2007. larXiv:0709.0497l [hep-ph]; 
V. Pilipp, Nucl. Phys. B 794 (2008) 154 [ arXiv: 0709 .32141 [hep-ph]]. 

[4] M. Beneke and S. Jager, Nucl. Phys. B 768 (2007) 51 |ar Xiv:hep-ph/0610322| ; 
A. Jain, I. Z. Rothstein and I. W. Stewart, arXiv: 0706 .33991 [hep-ph]. 

[5] G. Bell, Nucl. Phys. B 795 (2008) 1 |arXiv:0705.3T27l [hep-ph]]. 



28 



[6] G. Bell, PhD thesis, LMU Miinchen, 2006. larXiv:0705.3133l [hep-ph]. 

[7] K. G. Chetyrkin, M. Misiak and M. Munz, Nucl. Phys. B 520 (1998) 279 
[arXiv:hep-ph/9711280| . 



M. Gorbahn and U. Haisch, Nucl. Phys. B 713 (2005) 291 | arXiv:hep-ph/04110711 . 

[9] A. J. Buras, M. Gorbahn, U. Haisch and U. Nierste, JHEP 0611 (2006) 002 
[arXiv:hep-ph/06030~79l . 

[10] M. Gorbahn, private communication. 

[11] F. V. Tkachov, Phys. Lett. B 100 (1981) 65; 

K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B 192 (1981) 159. 



[12] R. Bonciani and A. Ferroglia, JHEP 0811 (2008) 065 [arXiv:0809.4687 [hep-ph]]; 
H. M. Asatrian, C. Greub and B. D. Pecjak, Phys. Rev. D 78 (2008) 114028 
jarXiv:0810.0987l [hep-ph]]. 

[13] M. Beneke, T. Huber and X. Q. Li, Nucl. Phys. B 811 (2009) 77 [ arXiv:0810.1230l 
[hep-ph]]. 

[14] E. Remiddi and J. A. M. Vermaseren, Int. J. Mod. Phys. A 15 (2000) 725 
[arXiv:hep-ph/99052"37| . 

[15] G. Bell, Nucl. Phys. B 812 (2009) 264 |arXiv:0810.5695l [hep-ph]]. 

[16] T. Huber. lirXiv:0901.2133 [hep-ph]. 

[17] P. Gambino, M. Gorbahn and U. Haisch, Nucl. Phys. B 673 (2003) 238 
[arXiv:hep-ph/03060'79| . 

[18] A. V. Efremov and A. V. Radyushkin, Phys. Lett. B 94 (1980) 245; 
G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22 (1980) 2157. 

[19] T. Becher, M. Neubert and B. D. Pecjak, Nucl. Phys. B 619 (2001) 538 
[arXiv:hep-ph/01022"l9| ; 

C. N. Burrell and A. R. Williamson, Phys. Rev. D 73 (2006) 114004 
|arXiv:hep-ph/0504024| . 

[20] R. J. Hill, T. Becher, S. J. Lee and M. Neubert, JHEP 0407 (2004) 081 
|arXiv:hep-ph/04042"l7| ; 



M. Beneke and D. Yang, Nucl. Phys. B 736 (2006) 34 |arXiv:hep-p h/0508250|. 
[21] S. J. Lee and M. Neubert, Phys. Rev. D 72 (2005) 094028 |arXiv:hep-ph/0509350] . 
[22] T. Becher and R. J. Hill, JHEP 0410 (2004) 055 [arXiv:hep-ph/0408344] ; 



G. G. Kirilin, arXiv:hep-ph/0508235, 



29 



[23] C. Bobeth, M. Misiak and J. Urban, Nucl. Phys. B 574 (2000) 291 
|arXiv:hep-ph/9910220| . 

[24] O. V. Tarasov, A. A. Vladimirov and A. Y. Zharkov, Phys. Lett. B 93 (1980) 429; 
S. A. Larin and J. A. M. Vermaseren, Phys. Lett. B 303 (1993) 334 
|arXiv:hep-ph/93022"08| . 

[25] M. Misiak and M. Munz, Phys. Lett. B 344 (1995) 308 | arXiv:hep-ph/9409454] ; 
K. G. Chetyrkin, M. Misiak and M. Munz, Nucl. Phys. B 518 (1998) 473 
|arXiv:hep-ph/9711266| . 

[26] E. Egorian and O. V. Tarasov, Teor. Mat. Fiz. 41 (1979) 26 [Theor. Math. Phys. 41 
(1979) 863]. 

[27] R. Tarrach, Nucl. Phys. B 183 (1981) 384; 

D. J. Broadhurst and A. G. Grozin, Phys. Rev. D 52 (1995) 4082 
[arXiv:hep-ph/9410240| . 



30 



